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We give a short introduction to the Analytic Perturbation Theory (APT) [T] 
and its generalization to Fractional powers — FAPT [51 [3] . We describe 
how to treat heavy-quark thresholds in FAPT and then show how to resum 
perturbative series in both the one-loop APT and FAPT. As an application 
we consider FAPT description of the Higgs boson decay H° — > bb. 



1. APT and FAPT in QCD 

In the standard QCD Perturbation Theory (PT) we know the Renor- 
malization Group (RG) equation da s [L]/dL = —a 2 — . . . for the effective 
coupling a s (Q 2 ) = a a [L)/(3 f with L = HQ 2 /A 2 ), (3 f = b (N f )/(4ir) = 
(11 — 2N f / 3) / (Airj^. Then the one- loop solution generates Landau pole 
singularity, a s [L] = 1/L. 

In the Analytic Perturbation Theory (APT) we have different effective 
couplings in Minkowskian (Radyushkin [3], and Krasnikov and Pivovarov 
[5]) and Euclidean (Shirkov and Solovtsov [1]) regions. In Euclidean domain, 
—q 2 = Q 2 , L = lnQ 2 /A 2 , APT generates the following set of images for the 
effective coupling and its n-th powers, {-4n[-^]} ne N' whereas in Minkowskian 
domain, q 2 = s, L s = lns/A 2 , it generates another set, {2l n [L s ]} neN . APT 
is based on the RG and causality that guaranties standard perturbative UV 
asymptotics and spectral properties. Power series Ylm^rnO^[L\ transforms 
into non-power series ^2 m d m A m [L] in APT. 

* Presented at the International Meeting "Excited QCD", February 8-14, 2009, Za- 
kopane (Poland) 

1 We use notations f(Q 2 ) and f[L] in order to specify the arguments we mean — 
squared momentum Q 2 or its logarithm L = ln(Q 2 /A 2 ), that is f[L] = /(A 2 ■ e L ) and 
A 2 is usually referred to Nf = 3 region. 
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By the analytization in APT for an observable f(Q 2 ) we mean the 
"Kallen-Lehman" representation 

Then in the one-loop approximation for the running coupling its spectral 
density is pi(cx) = 1/ ' \j ' L 2 + ir 2 and 



ML]= I ^ 2 dv= 1 I -— i , (2a) 



?l,(Z,l - / '^-da = larccos^^' ( 2b ) 
cr vr V 71 " + L s 



whereas analytic images of the higher powers (n > 2, n E N) are: 
A„[L]\ 1 / d\ n -Vii[L] 



2l n [L s y (n- 1)! V 



(3) 



In the standard QCD PT we have also: 

(i) the factorization procedure in QCD that gives rise to the appearance of 
logarithmic factors of the type: a v s [L] L; 

(ii) the RG evolution that generates evolution factors of the type: B(Q 2 ) = 
[Z(Q 2 )/Z(fj, 2 )] B(n 2 ), which reduce in the one-loop approximation to 
Z(Q 2 ) ~ a u s [L] with v = 70/(260) being a fractional number. 

All these means we need to construct analytic images of new functions: 
a s , a v s L m , .... 

In the one-loop approximation using recursive relation (|3|) we can obtain 
explicit expressions for ^[L] and %l u [L]: 

A m- 1 F{e~ L ,l-v) [n _ sin [^-^^(t^.j ,„ 
AAM -L» T(u) '^ [LJ - ^-l)( vr 2 + L 2 ) (-D/ 2 - W 

Here F(z, v) is the reduced Lerch transcendental function, which is an an- 
alytic function in v. They have very interesting properties, which we dis- 
cussed extensively in our previous papers [21 [7]. 

Construction of FAPT with fixed number of quark flavors, Nf, is a two- 
step procedure: we start with the perturbative result [a s ((5 2 )] iy , generate 
the spectral density p v {cr) using Eq. ([T]), and then obtain analytic couplings 
A U [L] and 2lj,[X] via Eqs. ([2]). Here Nf is fixed and factorized out. We can 



2 First indication that a special "analytization" procedure is needed to handle these 
logarithmic terms appeared in [6]. 
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proceed in the same manner for iVy-dependent quantities: [a s (Q 2 ; iVy)]" =>• 
p u (a; N f ) = p u [L a ; N f ] = p v {p)l(5) A U [L; N f ] and % v [L; N f ] — here N f 
is fixed, but not factorized out0 

Global version of FAPT, which takes into account heavy-quark thresh- 
olds, is constructed along the same lines but starting from global perturba- 
tive coupling [af lob (Q 2 )] , being a continuous function of Q 2 due to choosing 
different values of QCD scales Ay, corresponding to different values of Nf. 
We illustrate here the case of only one heavy-quark threshold at s = m 4 , 
corresponding to the transition Nf = 3 — > Nf = 4. Then we obtain the 
discontinuous spectral density 

pt » = 8(L a < L 4 ) p n [L a ; 3] + 6 (L 4 < L a ) p n [L a + A 4 ; 4] , (5) 

with L a = In (cr/Af), L f = In fmJ/A§) and A/ = In (A|/Aj) for / = 4, 

which is expressed in terms of fixed-flavor spectral densities with 3 and 
4 flavors, p n [L;3] and p n [L + A 4 ;4]. However it generates the continuous 
Minkowskian coupling 

2lf b [L s ] = (L s <L 4 ) (2t,[L s ; 3] - ^[L 4 ; 3] + %[L 4 + A 4 ; 4]) 

+ 0(L 4 <L s ) 2UL S + A 4 ;4]. (6) 

and the analytic Euclidean coupling ^ lob [L] (for more detail see in [7]). 



2. Resummation in the one-loop APT and FAPT 

We consider now the perturbative expansion of a typical physical quan- 
tity, like the Adler function and the ratio R, in the one-loop APT. Due 
to the limited space of our presentation we provide all formulas only for 
quantities in Minkowski region: 

oo 

n[L]=d + Y,d n ^n[L}- (7) 

n=l 

We suggest that there exist the generating function P(t) for coefficients 
d n = dn/dr. 

f'OQ f'OQ 

d n = P(t)t n ~ x dt with / P{t)dt = l. (8) 
Jo Jo 



3 Remind here that j3 s = b (N f )/(4ir). 
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To shorten our formulae, we use for the integral f °°f(t)P(t)dt the following 
notation: ((f(t)))p(t)- Then coefficients d n = d\ ((t n ~ 1 )} pu\ and as has 
been shown in [8] we have the exact result for the sum in (|7|) 



K[L]=d Q + d 1 ((% 1 [L-t]))p (t) 



(9) 



The integral in variable t here has a rigorous meaning, ensured by the finite- 
ness of the coupling 2li[t] < 1 and fast fall-off of the generating function 
P(t). 

In our previous publications [3 |9] we have constructed generalizations 
of these results, first, to the case of the global APT, when heavy-quark 
thresholds are taken into account. Then one starts with the series of the 
type (|7j) , where 2l„ [L] are substituted by their global analogs 2l| lob [L] (note 
that due to different normalizations of global couplings , 2l glob [L] ~ 2l rt [L] / (5f , 
the coefficients d n should be also changed). Then 



7£ glob [L] = d +d 1 ((9(L<L 4 ) A 4 2li[t] +21] 
+d 1 {{e(L>L 4 )% 



L--;3 

P3 



)>P(t) 



L+X 4 - — ;4 

Pa ■ 



))p(t) ; 



(10) 



where A 4 2l^[t] = % v L 4 + A 4 - t/fa; 4 



21, 



L 3 -t/p 3 ;3 



The second generalization has been obtained for the case of the global 
FAPT. Then the starting point is the series of the type Yl^=o dn ^n+A-L] an d 
the result of summation is a complete analog of Eq. (fl"0j) with substitutions 



P{t) =► P v {t) 



P 



1-x 



V X 



v-l 



dx 



1 — x 



(11) 



d => do%u[L], %i[L ~ t) => %i+u[L - t], and A 4 2l a [i] => A 4 2ti + „[i]. All 
needed formulas have been also obtained in parallel for the Euclidean case. 



3. Applications to Higgs boson decay 

Here we analyze the Higgs boson decay to a bb pair. For its width we 
have 



T(H -» bb) 



4V^7T 



M H R S {M 2 H ) with R S (M 2 H ) = m 2 b (M 2 H )R s (M 2 H ) (12) 



Resummation in FAPT: How many loops to calculate? 



■5 



and R s (s) is the i?-ratio for the scalar correlator, see for details in [2| I10j. 
In the one-loop FAPT this generates the following non-power expansion^: 

K 3 [L] = 3^, lfXSt[L] +d? $%t+» lL]\ , (13) 

I n>l J 

where fw?n = 8.45 GeV 2 is the RG-invariant of the one-loop m 2 (// 2 ) evo- 
lution ml(Q 2 ) = rh 2 {l) <°(Q 2 ) with u = 2j /b (5) = 1.04 and 70 is the 
quark-mass anomalous dimension (for a discussion — see in 




Fig. 1. Left: The relative errors A S N [L], N = 2,3 and 4, of the truncated FAPT in 

comparison with the exact summation result, Eq. (fl~5|) . Right: The width 

as a function of the Higgs boson mass Mjj in the resummed FAPT (solid line) . 



We take for the generating function P(t) the Lipatov-like model of [9] 
with {c = 2.4, = -0.52} 



l + (3 



Fs{t} - c(l + (3) e 



(14) 



It gives a very good prediction for d^ with n = 2, 3, 4, calculated in the QCD 
PT [10]: 7.50, 61.1, and 625 in comparison with 7.42, 62.3, and 620. Then 
we apply FAPT resummation technique to estimate how good is FAPT in 
approximating the whole sum 7£ S [L] in the range L £ [11.5,13.7] which 
corresponds to the range M H £ [60, 180] GeV 2 with Aq£= 3 = 189 MeV and 
2lf ob (m|) = 0.122. In this range we have (L 6 = ln(m 2 /A|)) 



ML] 
3m 2 1} 



2lfI° b [L] + ^ ((% 1+vo 



7T 



L+X 5 — 



+ A 6 2l 



» P s (15) 



4 Appearance of denominators ir n in association with the coefficients d n is due to d n 
normalization. 
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with P„ (t) defined via Eqs. (I14p and (llip . Now we analyze the accuracy of 
the truncated FAPT expressions 



K S [L;N] 



n - 2 

3m (1) 



N 



7T" 



■aCj£] 



(16) 



and compare them with the total sum 1Z S [L] in Eq. fjl5f) using relative errors 
A^[L] = 1 — 1Z S [L; N]/TZ S [L]. In the left panel of Fig. CD we show these errors 
for N = 2, N = 3, and iV = 4 in the analyzed range of L G [11, 13.8]. We see 
that already 1Z S [L;2] gives accuracy of the order of 2.5%, whereas 7£ S [L;3] 
of the order of 1%. That means that there is no need to calculate further 
corrections: at the level of accuracy of 1% it is quite enough to take into 
account only coefficients up to d%. This conclusion is stable with respect 
to the variation of parameters of the model P s (t) and is in a complete 
agreement with Kataev-Kim conclusion |llj . 



4. Conclusions 

In this report we described the resummation approach in the global 
versions of the one-loop APT and FAPT and argued that it produces finite 
answers, provided the generating function P(t) of perturbative coefficients 
d n is known. The main conclusion is: To achieve an accuracy of the order 
of 1% we do not need to calculate more than four loops and coefficients 
are needed only to estimate corresponding generating functions P(t). 
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